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Abstract 

The discrete Painleve I equation (dPi) is an integrable difference equation 
which has the classical first Painleve equation (Pi) as a continuum limit. dPi 
is believed to be integrable because it is the discrete isomonodromy condi- 
tion for an associated (single- valued) linear problem. In this paper, we derive 
higher-order difference equations as isomonodromy conditions that are asso- 
ciated to the same linear deformation problem. These form a hierarchy that 
may be compared to hierarchies of integrable ordinary differential equations 
(ODEs). We strengthen this comparison by continuum limit calculations 
that lead to equations in the Pi hierarchy. We propose that our difference 
equations are discrete versions of higher-order Painleve equations. 



1 Introduction 

Our aim is to derive higher-order versions of the equation 

an + /3 

Vn+i + Vn + Vn-i = 1- 7- I 1 ) 

Vn 

Equation ([l|) is known as the first discrete Painleve equation or dPi, because 
the scaled continuum limit y n — 1 + h 2 u(x), x = nh, with a = ri/i 5 , (3 = 



—3 + r 2 h 2 , 7 = 3-/3 — r 3 /i 4 , where the r« (i = 1, 2, 3) are constants, yields a 
scaled and translated version of the classical first Painleve equation (Pi) 



u" = Qu + x, 

as h — > 0. 

Pi is the simplest of six well known nonlinear second-order ODEs in the 
complex plane called the Painleve equations. Their characteristic property 
that all movable singularities of all solutions are poles is called the Painleve 
property. Painleve Q, Gambier 0, and Fuchs || identified them (under 
some mild conditions) as the only such equations with the Painleve property 
whose general solutions are new transcendental functions. 

Pi is believed to be integrable because it is the isomonodromy condition for 
an associated (single- valued) linear system of differential equations || . dPi is 
also an isomonodromy condition. Moreover, it possesses a discrete version of 
the Painleve property called the singularity confinement property proposed 
by Grammaticos, Ramani et al^ 

Joshi et al J7j derived the linear problem associated with dPi by starting 
with the Ablowitz-Ladik ||], [| scattering problem. The latter authors in 
turn based their scattering problem on the general AKNS problem given by 
Ablowitz, Kaup, Newell, and Segur [|l(| . 



It is well known that a single linear problem can give rise to a hierarchy 



of integrable ODEs ||11|| . In this paper, we announce an extension of such 
hierarchies to the discrete realm by deriving hierarchies of integrable difference 
equations. 

The plan of the paper is as follows. We recall the isomonodromy prob- 
lem for dPi in section 2 for completeness and then show how higher-order 
difference equations arise from it, giving examples up to the sixth-order. 

In section 3, we give continuum limits of the second and fourth-order dPi 
equations. In section 4, we examine a linear problem associated with Pi and 
show that the continuum limit found for the fourth-order dPi is the next 
equation in the Pi hierarchy. We propose that our hierarchy of difference 
equations is an integrable discrete version of such a hierarchy. 



2 Construction of the dPj hierarchy 

The linear problem associated with dPi |7j is 

a n w n+ i = \w n -w n -i, (2) 

-qT- = a n w n+ i + b n w n , (3) 

where w, a, a, b depend on a discrete variable n and w, a, b also depend on the 
continuous variable A (a, b rational in A). 



The aim of this section is to show that dPi is one of a whole range of 
higher order compatibility conditions arising from this system. 
The compatibility conditions of equations (|2|) and are 

b n +i - K-\ + A(- 5± ^ - — ) = 0, 

A 2 ,a ra +i fl n n A 1 a n+1 

( ) H (On+1 - »n) (a« On-1 + 1) = 0. 

Define 

a 

Pn = — and q n = b n + 6 n _i. 

We concentrate on results for p n . Corresponding results for q n can be obtained 
from the first compatibility condition above. 

The two compatibility conditions collapse down into a single condition in 
terms of p n 

Un+iPn+2 + {a n - A 2 )p„+i + (A 2 - a n )p n - a n _ip„_i + 2 = 0. (4) 

Take p n to be a finite Laurent expansion in A 

I 

Pn = X! -^\"A ! 
fc=0 

and substitute this expansion into (B). (Wherever convenient below we refer 
to Pjfc.n as Pfc.) We then find the following set of simultaneous equations for 
P P P 

1 0,n i 1 l,n i ■ l,n 

Pl,n ~ Pl,n+1 — 0, 
Pl-l,n ~ Pl-l,n+l — 0j 

ttn+l-Pfc,n+2 + «n-Pfc,n+l — «nPfc,n ~ &n-\Pk,n-\ + Pk-2,n — Pfc-2,n+l = 

for 2 < k < I, and 

QJn+l-f l,n+2 + a n-Pl,n+l — Oi n P\ jTl — Q! n -l-fl,n-l = 
ttn+l-Po,n+2 + «n-Po,n+l — ttn-Po,n — ttn-lPo.n-l + 2 = 

for fc = 0, 1. 

Note that the equation for Po is inhomogeneous (i.e. it has a nonzero 
forcing term) whereas the equation for Pi is homogeneous. Also note that 
the equations defining P*. for even and odd k may be separated into otherwise 
identical equations. These lead to two inconsistent compatibility conditions 
for a unless one of these two subsequences of P/. vanishes identically. In 



fact, Pi = for odd I is the only possibility due to the inhomogeneity of the 
equation for P . 

We are therefore left with the following system 

P2m,n ~ p2m,n+l = 0) 
&n+lP2k,n+2 + & n P 2 k, n +l ~ Ot n P2k,n ~ Ot n -iP 2 k,n-l + P2k-2,n ~ P2k-2,n+l = 

for 1 < k < m, and 

«n+l-Po,n+2 + «n-Po,n+l — «n-P(),n — &n-lPo,n-l + 2 = 

for k = 0, where we have taken / = 2m (m £ N). 

Below we list the solutions for different m up to m = 3. The q's, i = 
1, 2, 5 are constants. 

• m = yields the trivial linear nonautonomous equation 

a„c 2 = ci + Co(-l) n - n 

with 

P ,n = C 2 ^ 0. 

• m = 1 yields the second-order equation 

a n c 3 (a n -i + a n + a n+1 ) + a n c 2 = c\ + c (-l) n - n (5) 

with 

Po, n = c 2 + c 3 (a n + a„_i), 
P2,n = c 3 ± 0. 

In this case, we have recovered the more general version of dPi f|. 

• m = 2 yields the fourth-order equation 



a n c 4 (a n+ ia n+2 + a 2 n+1 + 2a n a n+ i + a ri „ia n _ 2 + o? n _ x + 2ct„a n _i 

+ a n -ia n+ i) + a„c 3 (a n _i + «„ + a n +i) + a n c 2 
= Cl + Co (-l) n -n (6) 

with 

Po,n = c 2 + c 3 (a n + a n -i) 

+c 4 (2a n a„,_i + a^_! + a n _ia n _ 2 + a n a n+1 + a 2 n ), 
P%n = c 3 + c 4 (a n + a n _i), 

P 4 ,n = C 4 ^ 0. 



• m = 3 yields the sixth-order equation 

a n C5(a n+ ia n+2 a n+ 3 + 2a n a n+1 a n+2 + a n -ia n+ ia n+2 + a n+1 
+a n+ ia 2 n+2 + 2a n+2 a 2 l+l + 3a n+ ia^ + 3a n a^ +1 + a n -ia n - 2 a n - 3 
+2a n a n _ia n _ 2 + a n _ 2 a n -ict„ + i + a? n _ x + a^_ 2 a„_i + 2a 2 n _ 1 a n _ 2 
+3a n _i«^ + 3a£_ia„ + 4a n _ia n a n+ i + a£ + a n _ia£ +1 
+«n-i a n+i) + «nC4(an+ia n .+2 + ai 2 +1 + 2a„,a: n+ i + a„_io; n _2 
+a^,_i + 2a n a n _i + a 2 n + a„-ia n+ i) + a„c 3 (a n _i + a n + a n +i) 
+a n c 2 = ci + c (-l) n - n 

with 

fo,n = c 2 + c 3 (a n + a n -i) + c 4 (2a n a n _i + a£_i + a n _ict„_ 2 

+a„a n+ i + a 2 ) + c 5 (a n a n+1 a n+2 + 2a n _ia„a; n+ i + 3a 2 a n -i 

2 2 2 3 3 2 

+3a n _ 1 a„ + a n _ 2 a„_i + 2a n+ \a n + a n + a n _ 1 + a n a n+1 
+a n _ia n _ 2 a n _ 3 + 2a n a n _ia n _ 2 + 2a^_ 1 a„_ 2 ), 

P2, n = c 3 + c 4 (a n + a n -i) 

+c 5 (2a„a n _i + a 2 n _ x + a„_ia n _ 2 + a n a n+1 + a 2 n ), 

P 4 , n = c 4 + c 5 (a n + a n _i), 

-^6,?! — C5 7^ 0- 

This process of solving the system for increasing m continues indefinitely. 
At each stage, the order of the compatibility condition increases by 2. 



3 Continuum limits 

In this section, we derive continuum 
the previous section. 



limits of difference equations found in 



3.1 The case m = 1 

Consider equation @. If Co = 0, we recover the equation often referred to 
as dPi, with Pi as one of its continuum limits. If Cq 7^ however, the term 
(—1)™ suggests an odd-even dependence in a n . This dependence must be 
taken into account to obtain a meaningful continuum limit. This leads to the 
transformation 

«2fe-i = u k , a 2k = v k . 



(This is similar to the limit pointed out by Grammaticos et al [12) for a 



discrete version of the second Painleve equation generalised with an additional 
(—1)" term.) 



In our search for a continuum limit, we use the substitutions 

2 

Uf. = 1 + hy(kh), Vk = z(kh), t = kh, C3 = h~ 3 , 

n 

and for ease of notation rename 

ci + c ci - Co + 1 c 2 2 
/i = , i/ = , cr = and p = . 

c 3 c 3 c 3 c 3 

Then we find 

a + u — 1 u + 1 1 2 z/ 1 2 

* = g 2~~ ^ + + ^ ~ 2 y * ~ 2^ 

+ ^ - ntj/ - z/y 3 )/i 3 + 0(/* 4 ). 



Using the scalings 

a 2 v 2 v 3 

u = 1 1 h cr r 2 h , 

P 4 4 2 4 

a = 2 + y + ^ 2 , 
2 

1/ = 1 + -r 4 ft, 
in the limit /i^Owe are left with 

2 

- 2y 3 + 2r i y 2 - 2r x ty + 2r 3 y + -nrjt + r x + 2r 2 = 0. 
This is a scaled and translated version of the second Painleve equation (Pn). 

3.2 The case m = 2 

For simplicity, we restrict ourselves to the case Co = 0. First, we illustrate the 
fact that second-order continuum limits are possible even though difference 
equation (||) is fourth order. Under the substitution 

a n = 1 + h 2 u(nh), t = nh, 

and the scalings 

10 + ^-^ + ^ = ^, 20 + ^ + ^ = r 2 h 2 , - = r,h\ 
C4 C4 C4 C4 C4 C4 



equation (0) becomes 



(10 + -)uu + (10 - -)u 2 + r 2 u + n + r 3 t = 

C4 C4 



in the limit h —> 0. This is a scaled and translated version of Pi. However, 
this case restricts the four degrees of freedom contained in the parameters of 
equation (||) to three. 

Scalings that maintain the full four degrees of freedom, i.e. 

in . 3c 3 Ci c 2 6 3c 3 C\ 4 
10 H 1 —r x h\ 20 H 1 = r 2 h , 

C4 C4 C4 C4 C4 

10 + - = r 3 /i 2 , - = r 4 h 7 , 

C4 C4 

lead, in the limit h — > 0, to the fourth-order equation 

«tttt + 5(m 4 ) 2 + 10wu tt + r 3 u tt + 10u 3 + 3r 3 u 2 + r 2 u + n + r 4 t = 0. (7) 
We discuss the significance of this equation in the next section. 



4 Discussion 

The aim of this section is to show that equation (0), the continuum limit for 
the case m = 2, cq = 0, is the next equation in the Pi hierarchy. 

The continuum limit equations found for the case m = 1 are Pi and Pn. 
These equations are considered integrable as they are isomonodromy condi- 
tions for a single-valued linear problem. Furthermore, they may be viewed as 
equations lying at the base of two integrable hierarchies of equations. Con- 
sider the following isomonodromy problem associated with Pi 

-u(x)-e J) V ' Vk= (c -a) V ' (8) 
where u is the potential, k is the eigenvalue and 

B = (2r 2 + 4:u)k-8k 3 , 

C = -(2u xx + 2r 2 u + Au 2 )k + (-2r 2 + Au)k 3 + 8k 5 . 

Pi is a compatibility condition for system (§) with these choices of A, B, C. 
A hierarchy of compatibility conditions arises when A, B and C are expanded 
to higher degree in k. The next self-consistent expansions are 

A ( S ° , 3s Si S 3 

A = -{—u xxx + —uu x + —u x )k + —u x k, 

B = ^u xx + ^u 2 + S ^u-s 2 )k-C^u + Sl )k 3 + sok 5 , (9) 

o o 2 Z 

C = {^uu xx - ^u 2 x + ^-u 3 + ^-u 2 -2x + s 3 )k 

+ {^-u xx + ~~u 2 + ~~u + s 2 )k 3 + {-^-u + si)k 5 - s k 7 , 

o o 2 2 



where sq,si, sg, S3 are constants. A straight forward calculation shows that 
equation (|7|) is a compatibility condition of (§) with A, .B, C expanded as in 
(D (renaming s = -g, *i = s 2 = s 3 = 
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